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Random walk on discrete lattices is a fundamental model in physics that forms the basis for 
our understanding of transport and diffusion process. In this work, we study unreachability in 
networks, i. e, the number of nodes not visited by any walkers until some Hnite time. We show that 
for the case of multiple random walkers on scale-free networks, the fraction of sites not visited is 
well approximated by a stretched exponential function. We also discuss some preliminary results 
for distinct sites visited on time-varying networks. 
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I. INTRODUCTION 

Random walk was introduced more than 100 years ago 
by Karl Pearson [I| and over the decades this idea has 
formed the basis for our understandingof transport and 
diffusion processes in physical systems [ 4 I . Random walk 
model continues to be relevant today in many applica¬ 
tions ranging from physics and computer sciences to biol¬ 
ogy i- In particular, random walks on disordered lattice 
such as complex networks Q form the basis for modern 
applications such as web search engines @ and recom- 
mender systems 0. In addition, search and discovery on 
complex networks [7| , information spreading on networks 
0 , clustering algorithms @ and extreme events on net¬ 
works [ 1 ^ are other areas in which random walks play a 
crucial role. 

One class of relevant application is in the design of 
algorithms for efficient search and discovery process on 
complex networks. What route should a packet of in¬ 
formation on information networks travel from origin to 
its destination ? How to quickly find new webpages on 
the world wide web ? All these questions are similar in 
spirit to the general problem of search for information or 
other resources on complex networks. The process of dis¬ 
covery of resources through random walking on networks 
is useful as a benchmark to against which to measure 
the efficiency of other search and discover algrithms. In 
many practical networked systems several variants of ran¬ 
dom walks are used as algorithm to explore a network, 
TCP/IP being a specific example m- 

As a statistical physics problem, random walks on dis¬ 
crete lattices has been studied for several decades espe¬ 
cially to understand first passage time distribution |12l |. 
distance covered by walkers , extreme events M , distinct 
sites visited and site occupancy, flux-fluctuation relations 
0 and node centrality in complex networks Most 
of these measures have direct bearing on the questions 
related to efficient search and information diffusion on 
networks. Many results for random walks on regular d- 
dimensional lattices are known [l^ . such as for instance 
mean number of distinct sites visited by W random walk¬ 
ers [I- On regular d-dimensional lattices and for short 
times, the number of distinct sites visited f{t) oc and 


asymptotically behaves f{t) oc Wt for d > 3. In general, 
Ref. [13 identifies three distinct temporal regimes with 
different behaviours for fr{t). 

On the contrary, if the underlying topology is not a 
regular lattice but instead forms a complex network, not 
many results are known for distinct sites visited by mul¬ 
tiple random walkers. Distinct visited sites are those vis¬ 
ited by at least one walker until time t. For the case of a 
single walker on scale-free networks, it has been numer¬ 
ically shown that the fraction of distinct sites visited in 
t time steps is fr{t) = t and as t —>■ co, /r(t) —>■ 1 due 
to finite size of network [l^ . In a small world network, 
fr-{t) displays a cross-over from fr{t) = to = t 

behaviour depending on whether the walker has managed 
to hit a short-cut in the SW network or not 0. Another 
related point of view is to consider the probability for the 
paths taken by the random walkers on networks and it 
has been shown that only three distinct distributions are 
possible [Tsil . 

Recently, exact distribution for the number of distinct 
sites visited by many random walkers in a one dimen¬ 
sional lattice has been obtained [l^. What happens in 
the case of multiple random walkers on a network ? In 
this work, we study unreachability furit), i.e, the frac¬ 
tion of sites not visited by any of the W random walkers 
on scale-free network until time t. We show that this 
fraction fur(t) can be modelled as a stretched exponen¬ 
tial function. If A/, is the number of nodes reached in 
time t, then fur{t) is related to fraction of distinct nodes 
fr{t) = Nr/N visited by the random walkers in time t 
through the relation fur{t) = 1 — fr{t). 

Further, we discuss the results obtained in the case 
when the networks are rewired. It is known that many 
real-life networks are time-varying [l^ and if the time 
scales over which network varies is comparable to the 
time scale of dynamics taking place on them then it be¬ 
comes necessary to account for temporal variation of net¬ 
works as well. We consider a case in which some specific 
percentage of nodes are rewired and we show that net¬ 
work reachability is enhanced by the process of rewiring 
the edges. 
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FIG. 1. Unreachability of random walkers on scale-free net¬ 
works plotted as a function of a; = WtjN. All the dots were 
obtained through direct random walk simulations. The solid 
(black) curve is the best fit line. Note the excellent scaling 
collapse. 

II. RANDOM WALKS ON NETWORKS 

Random walks on networks is a straightforward gen¬ 
eralisation of one-dimensional lattice random walk to a 
network setting. Thus, a walker on f-th node with de¬ 
gree di at discrete time n can choose to hop to any of 
its neighbours with the probability l/di. In this work, 
we consider a scale-free network with N nodes and E 
edges on which W independent random walkers are exe¬ 
cuting random walk. This dynamics can be described by 
a master equation of the form, 

P,,(n + 1) = ^^ P.fc(n) (1) 

k 

where Pij{n) is the occupation probability for a walker 
starting *-th node at time n = 0 and reaching j-th node 
at time n-|-l. In this, Aij represents the elements of a ad¬ 
jacency matrix which carries the information about con¬ 
nectivity of the network, Aij = 1 if i-th and j-th nodes 
are connected and 0 otherwise. Starting from this equa¬ 
tion, many statistical quantities related to the random 
walk such as the first passage time can be calculated. 

In Fig. [1] we show the fraction of sites not reached by 
random walkers in scale-free networks with degree distri¬ 
bution P{k) ^ where 7 ~ 2.2. At time t = 0, all the 
random walkers were placed on the same node designated 
as zeroth node. This figure shows the result for six cases 
with different network sizes N and number of walkers W. 
For all the cases, we note that the scaled parameter of in¬ 
terest is cc = Wt/N and fur{t) = 1 — Nr/N. This choice 
of scaling the variables leads to a data collapse as shown 
in Fig. [H The unreachability is well approximated by a 
stretched exponential function of the form, 

fur(t) = e-^ (2) 



FIG. 2 . Same as in Figure [T] except that the it is plotted as a 
function of {Wt/N)^. In this semi-log plot, we obtain a good 
straight line. The black curve is the best fit line. 

where x = Wt/N and B r; 0.75 and /3 « 0.8 and these 
numerical values have been obtained through a fit to the 
data shown in Fig. [T] As cc << 1, we can expand the 
right hand side Eq. [5] to obtain fur{t) R! 1 — Bx^. Thus, 
the fraction of distinct sites visited is fr{t) R! x^. We note 
that in Ref. [l^ it was argued based on numerical results 
that for the case of one walker the distinct sites reached 
is given by fr{t) = t/N for t << N, implying that the 
parameter /3 = 1. Physically, the stretched exponential 
function conveys that the distinct number of sites visited 
will equal the number of sites in a network only at infi¬ 
nite time. This is also in agreement with the asymptotic 
result that the first passage time distribution for a ran¬ 
dom walker displays an exponentially decaying tail M- 
In Fig. m the same results as in Fig. [T] are shown, but 
now as a function of x^. As expected, in the semi-log 
scale, we obtain a linearly decaying curve confirming the 
scaling observed for unreachability. Thus, unreachabil¬ 
ity is a finite time effect in connected networks since as 
f —>■ 00 all nodes are eventually reached. 

III. SITES VISITED ON A REWIRED 
NETWORK 

In this section, we present preliminary results related 
to the number of distinct sites visited on a periodically 
rewired network. This problem is relevant in the case of 
search on networks that are subject to rewiring of the 
edges. Consider the case of web graph whose nodes (web 
pages) get rewired almost all the time. In practice, it is 
useful to consider snapshots of a network like web graph 
at periodic intervals to monitor how random walk proper¬ 
ties are modified by the changes in the network structure. 

In the spirit of this argument, we consider the dy¬ 
namics of W random walkers on a scale-free network 
and the rewiring algorithm proceeds as follows. Con¬ 
sider a fully connected scale-free network with N nodes. 
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FIG. 3. Reachability on a rewired network. Results of ran¬ 
dom walk simulation on a network with 2000 nodes, W = 20 
walkers is shown. Results for various rewiring periods are 
displayed. 

We rewire p percentage of its nodes after every T time 
steps. For every t = nT, where n > 0 is an integer, 
we randomly choose two nodes, with probabilities gi{k) 
and 32 (fc) respectively, and exchange their edges. The 
probability of choosing a node with degree k is taken to 
be gi{k) = k°‘,i = 1,2. In this, ai = 0 denotes ran¬ 
dom choice of nodes without regard to their degree. For 
Ui > 0, hubs are preferentially chosen and for < 0 leaf 
nodes are preferentially chosen for rewiring. 

In Figure[3]we show the result for the fraction of nodes 
reached in t when the network is rewired after every T 
time steps. At each of the rewiring instants, p = 0.7% of 
the nodes were rewired. The value of ai = 1 and 02 = —3 
implies that every rewiring involves exchange of edges 
between one hub and one leaf node. As the figure reveals, 
more rewiring leads to more reachability of the network 
in comparison with the case of static network without 
rewiring. As T increases the time interval between two 
consecutive rewiring instances increases. Hence, as we 
would expect, reachability is less and tends to merge with 
the static network result. 

In general, increasing the value of rewiring percentage 
p also increases the reachability of the network in com¬ 
parison with that of the static network. Indeed, even 
if oi = 02 = 0 leading to random choice of nodes of 
rewiring, there is enhancement in the reachability of the 
network. This shown in Fig. |T] which reveals that peri¬ 
odic random rewiring also leads to enhanced reachability. 
However, if targetted rewiring is performed by selectively 
choosing hub and leaf, then reachability is even more en¬ 
hanced compared to the random rewiring case. 

Finally, we comment on the physical basis of the effect 
arising due to rewiring of the network. In a scale-free 


network there are few hubs but a large number of leaf 
nodes. Random walkers are most likely to be exploring 
the domain of large number of leaf hubs. However, when 
a rewiring is effected, two nodes are chosen and their 



FIG. 4. Reachability on a rewired network. Results of ran¬ 
dom walk simulation on a network with 2000 nodes, IF = 20 
walkers is shown. Results for random rewiring is compared 
with the case for hub-leaf rewiring. 


edges are exchanged. If the nodes happen to be a hub 
and a leaf node, then random walkers on leaf node sud¬ 
denly find themselves on a hub and can quickly explore 
territory previously unavailable when they were on a leaf 
node. As more rewiring takes place this effect leads to 
enhanced reachability on the network. When nodes are 
randomly chosen the enhacement is still seen but not as 
much as in the case of hub-leaf rewiring. This is because 
in a scale-free network most nodes being leaf nodes, any 
exchange of edges between them does not lead to en¬ 
hanced reachability. For the same reason, rewiring is 
unlikely to lead to any enhancements in reachability for 
Erdos-Renyi networks . 

In summary, we have studied the question of reacha¬ 
bility of network nodes when multiple and independent 
random walkers are exploring a scale-free network. We 
show through numerical simulations that the unreacha¬ 
bility of the network, defined as the fraction of sites not 
visited until time t, can be represented by a stretched 
exponential function with x = Wt/N being the scaled 
parameter. As Fig. [1] shows, we obtain a good data 
collapse for various network sizes and numbers of walk¬ 
ers. Further, we also discuss some preliminary results 
on the effect of rewiring on the reachability of network 
nodes. We show that rewiring leads to enhancement in 
the number of distinct sites visited by random walkers. 
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